function [pval2s] = f_pval2S(tRet, T)

nvar   = length(tRet);
pval2s = nan(1,nvar);

if length(T)<nvar
    T = repmat(T,1,nvar);
end

for i=1:nvar
    pval2s(1,i) = 2*((tcdf(-abs(tRet(i)),T(i)-1))); % [H0: mu=0 H1: mu \= 0]
end